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Learning is but an adjunct to ourself 
And where we are our learning likewise is. 

- W. Shakespeare 



Abstract 



The paper studies ways in which the sets of a partition of a lattice in R" become regular model 
sets. The main theorem gives equivalent conditions which assure that a matrix substitution 
system on a lattice in R" gives rise to regular model sets (based on p-adic-like internal spaces) , 
and hence to pure point diffractive sets. The methods developed here are used to show that the 
n— dimensional chair tiling and the sphinx tiling are pure point diffractive. 

1 Introduction 

There have been two very successful approaches to building discrete mathematical structures with 
long-range aperiodic order. These are the substitution methods, notably symbolic substitutions and 
tiling substitutions, and the cut and project method. In the first case the structure is typically 
generated by successive substitution from a finite starting configuration. In the second it typically 
appears in one shot as the (partial) projection of a periodic structure in some "higher" dimensional 
embedding space. 

The principal focus in this paper is the relationship between matrix substitution systems on a 
lattice and a naturally related cut and project formalism. We start with a partition of a lattice L 
in K™ into a finite number of point sets U — (Ui, . . . , U m ) and a finite set of substitution rules $ 
which are affine inflations and under which U is invariant. The main theorem (Theorem |J) provides 
conditions on <f> which are equivalent to Ui, . . . , U m being regular model sets (i.e. cut and project 
sets) . One of the characterizations (modular coincidence) affords a simple computational approach to 
testing for model sets. In a later section we go beyond the context of substitution systems and provide 
an alternative characterization (Theorem |J) of model sets. We use both types of characterization in 
showing that the sphinx and n-dimensional chair tilings are based on model sets. 
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The connection between substitution systems and cut and project sets is nothing new, e.g. the 
Fibonacci chain is often described in terms of a cut by a strip through Z 2 , and the klotz construction 
of Kramer et al. || is a sophisticated elaboration of the same idea. Nonetheless, substitution systems 
and cut and project sets are not different formulations of the same thing, and the relationship between 
them remains inadequately understood. 

In the early study of aperiodic order, the cut and project formalism was always based on projection 
into R n from a lattice in some higher space R™ x R p , the projection being controlled by a compact set 
W ClM. p . However, it was already implicit in the much earlier work of Y. Meyer jl^] that R p can be 
replaced by any locally compact abelian group H and W C H by any compact set with non-empty 
interior, and the projection method still produces discrete aperiodic sets with diffractive properties 
(hence long-range order). Meyer's terminology for such sets was "model sets" and we use it here in 
deference to its priority and to emphasize the greater generality of the internal space H. Model sets 
have been studied in detail in (h], 13, [l4| [l5| [l(| |l7). The relevance of more general internal spaces 
to tiling theory and symbolic substitutions was made explicit in Q where p-adic and mixed p-adic 
and real spaces naturally appear. 

One of the important features of making the connection to model sets is that once it is established, 
pure point diffractivity is assured (see Theorem ^| for a precise statement of this). This type of 
information is generally quite difficult to obtain. For example, our results shown here prove that the 
ri-dimensional chair tiling and the 2-dimensional sphinx tiling are pure point diffractive. The former 
has been established for n = 2 previously [Q, [l8| . The latter is claimed in [|l8) as being provable by a 
geometric form of "coincidence" established there (see below for more on the concept of coincidence). 

The p-adic type internal spaces occur when the aperiodic set in question is based on the points 
of a lattice and its sublattices in R". An important class of examples of this type arises from the 
equal length symbolic substitution systems. Suppose that A = {a%, . . . , a m } is a finite alphabet with 
associated monoid of words A* , and we are given a primitive substitution a : A — > A* for which 
the length / of each of the words a(a,i) is the same. This substitution leads to a tiling of R of tiles 
of equal length, say equal to 1. Matching the coordinate of the left end of each tile with its tile type 
Oi, we obtain a partition U\ U . . . U U m of Z, and a may be viewed as comprised of a set of affine 
mappings x t— ► Ix + v where v G Z. 

A lot more is known about equal length substitutions than the arbitrary ones, a particularly 
important example of this being Dekking's criterion for diffraction ||. An equal length aperiodic 
tiling is pure point diffractive if and only if it admits a coincidence {a is said to admit a coincidence 
if there is a k, 1 < k < l n , for which the fcth letter of each word a n (ai) for some n is the same). 

In this paper we prove a related result, but this time the dimension is not restricted. Namely, 
there is a notion of coincidence (in fact there are two such notions) and either of these is equivalent 
to the sets U\, . . . , U m being regular model sets. One of the criteria for coincidence that we give is 
a straightforward algorithm and thus in principle is computable. 

As we have already pointed out, a consequence of our result is that coincidence implies the pure 
point diffractivity of Ui, . . . , U m . We do not know yet to what extent the condition is equivalent to 
pure point diffractivity. 

The setting of the paper is entirely at the level of point sets, so necessarily the strong conditions 
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implicit in the tiling situation are replaced here by a corresponding algebraic condition on the matrix 
substitution system: the Perron-Frobenius eigenvalue of the substitution system should equal its 
inflation constant. This is in fact a compatibility condition which is necessary for the model set 
connection to exist. This condition, not surprisingly, has occurred elsewhere in the literature (see for 
instance, Jll], ) . The important result that gets the process off the ground is Theorem [j], which is 
largely due to Martin Schlottmann. 

Matrix substitution systems, treated at the level of point sets, have recently appeared in Lagarias 
and Wang under the name of self-replicating Delone sets. In that paper, point sets X are not 
restricted to lattices and the principal question revolves around the interesting question of existence 
of tilings of W 1 by translations of certain prototiles for which the points of X are the appropriate 
translational vectors. Also related to our paper is the study of sets of affine mappings in the context 
of lattice tilings (see [[l9| for a nice recent survey on this). In relation to our paper, the situation there 
corresponds to the lxl matrix substitution systems and the problems become entirely different. 
Since the tilings there are lattice tilings, the whole issue of model sets and diffracion is trivial, and 
the issues lie more around the complex nature of the tiles themselves. 

2 Definitions and Notation 

Let X be a nonempty set. For m £ Z+, an to x to matrix function system (MFS) on X is an m x in 
matrix <!> = ($y ■), where each $^ is a set (possibly empty) of mappings X to X. 

The corresponding matrix S(<5>) :— (card^^))^ is called the substitution matrix of <I>. The MFS 
is primitive if S*($) is primitive, i.e. there is an I > for which S(<f>) 1 has no zero entries. 

In this paper we deal only with MFSs which are finite in the sense that card($ij) < oo for all 
Of particular importance are the Perron-Frobenius (PF) eigenvalue and the corresponding PF 
eigenvector (unique up to a scalar factor) of S(&). We will also have use for the incidence matrix 
/($) of <&, which is defined by 




if card($y) ^ , 
else . 



Let P(X) be the set of subsets of X. Any MFS induces a mapping on P(X) m by 

(1) 

which we call the substitution determined by We sometimes write ^ij(Uj) to mean f(Uj). 

In the sequel, X will be a lattice L in R™ and the mappings of $ will always be affine linear 
mappings of the form x i— > Qx + a, where Q £ End^(i) is the same for all the maps. Such maps 
extend to M. n . For any affine mapping / : x i— > Qx + b on L we denote the translational part, 6, of 
/ by t(f). We say that /, g £ $ are congruent mod QL if t(f ) = t(g) mod QL. This equivalence 
relation partitions $ into congruence classes. For a £ L, $[a] := { / £ (J i . | t(f) = a modQL}. 
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We say that $ admits a coincidence if there is an i, 1 < i < to, for which Hjli ^ $j i- e - the 
same map appears in every set of the i-th row for some i. Furthermore, if $ M [a] is contained entirely 
in one row of the MFS (<I> M ) for some M > 0, a G L, then we say that (t/,<J>) admits a modular 
coincidence. 

Let $,lbemxffl MFSs on X. Then we can compose them : 

*o$ = ((*o$) ii ), (2) 

where (* o = IX=i*«°*W and ^ fe o := | < 9 ° ^ ^ * J > Evidently 

5(* o $) < (see @ for the definition of the partial order). 

For an m x m MFS $, we say that U := \U\,--- , U m ] T <E P(X) m is a /ixed pomi of $ if $ [7 = U. 

3 Substitution Systems on Lattices 

Let L be a lattice in R™. A mapping Q £ Endz(L) is an inflation for L if det Q ^ and 

oo 

n ^ = w- ( 3 ) 

Let Q be an inflation. Then q := | dctQ| = \L : QL] > 1. We define the Q-adic completion 

L = L^ = lim L/Q k L (4) 

* — 

of L. L will be supplied with the usual topology of a profinite group. In particular, the cosets 
a + Q L, a G L, k = 0, 1, 2, • ■ • , form a basis of open sets of L and each of these cosets is both open 
and closed. When we use the word coset in this paper, we mean either a coset of the form a + Q k L 
in L or a + Q k L in L, according to the context. An important observation is that any two cosets in 
L are either disjoint or one is contained in the other. The same applies to cosets of L. 
We let fi denote Haar measure on L, normalized so that = 1. Thus for cosets, 

We also have need of the metric d on L defined via the standard norm: 

1 

g. 

From Plfclo Q k L = {0}, we conclude that the mapping x i— > {x mod Q k L}k embeds L in L. We 
identify L with its image in L. Note that L is the closure of L, whence the notation. 

An affine lattice substitution system on L with inflation Q is a pair (U, $) consisting of disjoint 
subsets {Ui}™ =1 of L and an m x to MFS $ on L for which U = [Ui, • • • , U m ] T is a fixed point of 
i.e. 

m 

Ui=\J (J /(C/,), i = l,...,m, (7) 



,,||:=J_ if x e Q fc L \ Q fc+1 L, 1 1 1 1 = . (6) 
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where the maps of $ are affine mappings of the form x i— » Qx + a, a G L , and in which the unions in 
(Q) are disjoint]] In this paper all our matrix substitution systems are composed of afhne mappings 
on a lattice and we often drop the words 'affine lattice', speaking simply of substitution systems. 

We say that the substitution system (U , $) is primitive if $ is primitive. A second substitution 
system (U' ,^f) is called equivalent to (U,&) if U' = U, $ and $ have the same inflation, and 
have the same PF-eigenvalue and right PF-eigenvector (up to scalar factor). 

Let (U, <J>) be a substitution system on L. Identifying L as a dense subgroup of L, we have a 
unique extension of <!> to a MFS on L in the obvious way. Thus if / € $y and / : x i— > + a, then 
this formula defines a mapping on L, to which we give the same name. Note that / is a contraction 
on L, since ||Qx|| = ^||a;|| for all x G L. Thus $ determines a multi-component iterated function 
system on L. Furthermore defining the compact subsets 

Wi~U~i, i=l,-..,m, (8) 

and using (|7|) and the continuity of the mapping, we have 

m 

W, = U U W). *' = !,••■,"», (9) 

which shows that W = [W\, ■ ■ ■ , W m ] T is the unique attractor of $ (see ||, §)• 

We call (W, $) the associated Q-adic system. We cannot expect in general that the decomposi- 
tion in (H) will be disjoint, so we do not call (W, $) a substitution system. 

For X,Y G M", we write 

X < Y" if X, < Y, for all 1 < i < n 

X < Y if Xi < Yi for all 1 < i < n . 

Similarly, for A,B e M„(K) 

A<B if Aij < Bij for all l<i,j<n , . 

A < B if Aij < Bjj for all 1 < i, j < n . ^ ^ 

We begin by recalling a couple of results from the Perron-Frobenius theory. 

Lemma 1 Let A be a non-negative primitive matrix with PF-eigenvalue A. IfO< XX < AX, then 
AX = XX. 

PROOF : We can assume X ^ 0. Since < XX and A > 0, X > 0. Let X' > be a PF right- 
eigenvector of A. Let a = max{ | 1 < i < m }. Then X < aX' and X is not strictly less than 
aX'. Claim X = aX'. If X ^ aV, then < A N {aX' - X) = aX N X' - A N X for some AT, 
since A is primitive. So X N X < A N X < aX N X', i.e. X < aX' . This is a contradiction. Therefore 
AX = XX. □ 



1 In the case that one has unions (|7|) which are not disjoint there arises the natural question of the mulitplicities 
of points, or more generally densities of points. For more on this see [j^, Elj. 
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Lemma 2 Let X be the PF- eigenvalue of the non-negative primitive matrix A and y, be an eigenvalue 
of a matrix B where < B < A. If A ^= B, then < A. 

PROOF : Let Y be a right eigenvector for eigenvalue fj, of B, with Y = [Yi, ■ ■ ■ ,Y m ) T . Let Y = 
[|Yi|, • • • , |y m |] T ^ 0. Then \fx\Y < BY < AY. Let X~ T be a positive left eigenvector for A with 
PF-eigenvalue A. So \fi\X T Y < if BY < if AY = XlfT. This shows that < A. If = A, 
then XY < AY. By Lemma [j], XY = AY. Since A is a primitive matrix, X m Y — A m Y > for some 
to. So Y > 0. From AF < BY < AY = XY, we have AY = BY. Therefore A = B. □ 

Lemma 3 Let (!/,$) be a primitive substitution system. Then for all I = 1,2, ([/,$') is a 
primitive substitution system. 

PROOF : Let i,j, k 6 {1,2, ■ ■ ■ ,m}. All the maps g £ have domain Uk and disjoint images in 
Ui. Moreover all the mappings g are injective. Likewise all the maps / of $jy have domain Uj and 
disjoint images in Uk- Thus all the maps g o f £ <& ik o have domain Uj and disjoint images in 
Ui. Furthermore $ 2 [7 = $(${/) = &{U) = U. So (U, $ 2 ) is a substitution system. The argument 
extends in the same way to (U, The statement on primitivity is clear. □ 

Theorem 1 Let (U,$>) be a primitive substitution system with inflation Q on L. Let (W, $) be 
the corresponding associated Q-adic system. Suppose that the PF-eigenvalue of S($) is detQ and 
L = \J™ 1 W i . Then 

(i) S(&) = (S(<P)y, r > 1; 

(ii) »(Wi) = ± YJj n =1 {S{^ r )) l3 y{W 3 ), for all i = l,---,m, r> 1; 
(Hi) For all i = !,-■■, m, Wi^ and fi(dWi) = 0. 

PROOF : For every measurable set E C L and all / £ <£>„ , fi(f(E)) = ji(Q(E) + a) = . d ^ t g, fi(E), 
where / : x i— > Qa; + a. In particular, /i(/(Wj)) = - Wj, where Wj :— (i(Wj) and q = | detQ|. We 
obtain 

^ - H "7 card (( $r )*i) W J 

from (|). 

Let w = [wi, ■ • ■ , w m ] T . Since [J™ -j Wi = L, the Baire category theorem assures us that for at 
least one i, 

Wi+ (11) 
and then the primitivity gives this for all i. So w > and 

w < —S(<P r )w < —S($Yw, for any r > 1 . (12) 

Since the PF-eigenvalue of S , ($) r is q r = | detQ| r and S , ($) r is primitive, we have from Lemma [j] 
that 

w = —S(<P r )w = —S($) r w, for any r > 1 . (13) 
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The positivity of w together with S($ r ) < S($) r shows that S($ r ) = S*($) r . This proves (i) and 
(ii). 

o . . 

Fix any i £ {1, • • •, m}, let Wi contain a basis open set a + Q r L with some r S Z>o by (|ll|). Since 

o 

(t/, $ r ) is a substitution system, a + Q r L C Wi C Wj = U^Li(* r )ij^"- In particular, (a + Q r L) n 
^(Wfc) ^ for some k £ {1, • ■ • , to} and some g G ($ r )ife. However g(L) = 6 + Q r L for some b 6 £, 
so (a + Q r L) n (6 + Q r Z) ^ 0. This means a + Q r L = b + Q r L. Thus 

g(W k ) C = a + Q r I . (14) 

o o 

For all / e (& r )ij, j € {1, 2, • • • , m}, / is clearly an open map, so U^i(^ r )^"(^j") c • Thus 

\J{® Uj(W s )J \Wi 



c 



3=1 

Note that due to ( pr) ) at least one 5 in ( < I >, ')ij does not contribute to the relation (|l5|). 

Let Vi :— fi(dWi), i = 1, • • • ,m and v :— [vi, ■ ■ ■ 7 v m ] T . So v < -^S(& r )v. Actually, by what we 
just said, 

< v < — S'v < — S($ r )v = — S($Yv, (16) 
q r q r q r 

where S' < S($) r ,S' ^ S(<Z>) r . Now applying the Lemma again we obtain equality throughout (|l(]) . 
But by Lemma |2] the eigenvalues of -^S' are strictly less in absolute value than the PF-eigenvalue 
of ^S , (<£>) r , which is 1. This forces v = 0, and hence /i(dWi) = 0, i = 1, ■ ■ ■ , m. □ 

In the sequel, the central concern is to relate the sets Ui and the sets A, :— Wi n L. Clearly 
Ai D Ui. The next lemma groups a circle of ideas that relate this question to the boundaries and 
interiors of the Wi. 

Lemma 4 Let Ui,i = 1, • • • , m, be point sets of the lattice L in R™ . Let Q be an inflation of L and 
identify L with its image in its Q-adic completion L. Define Wi :— Ui in L and Ai :— Wi D L. 



(i) IfUi,..., U m are disjoint and ^,(Ai\Ui) = for all i = 1, ■ • • ,m, then Wi fl Wj= for all 

o o 

(ii) Lf L = {J™ =1 Ui and Wi fl Wj= for all i ^ j, where i,j 6 {1, ■■■,m}, then Ai\Ui C 
U^i for all i — 1, • • ■ , to. 



(Hi) Ifn(dWj) = for all j = !,■■■, m and A l \U l C (T 1 dWj, then /i(A,\^) = 0. 
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PROOF: (i) Suppose there are i, j € {1, • • • , to} with W n Wjj^ 0. We can choose a <E (W t fl Wj) H L, 

o o o o 

since L is dense in L and Wi fl Wj is open. Choose k E Z + so that a + q k L dWi fl Wj. Note that 
a + q k L C Aj n Aj. Then 

m 

U( A *Wi) 2 ((a + g fe L)\;7 J )u((a + g fe L)\C/ 3 ) 

i=l 

= (a + g fe L) \ (C/i n Uj) 

= a + q k L, since the Ui, i = l,---,m, are disjoint . 

So 

m m 
i=l i=l 

> //(a + g fe L) 

> o, 

contrary to assumption. 

o o 

(ii) Assume Wi fl Wj — for all i ^ j. For any i G {1, • • • , m}, 

m 

(A 4 \C/ 4 ) C (U^)nWi, since L = |J (7, 

j^j i=l 
m 

C U(Wj n VK,) C |J since W 4 n W,- = for all i^j. 

(hi) Obvious. □ 



4 Model Sets 

Let us recall the notion of a model set (or cut and project set). A cut and project scheme (CPS) 
consists of a collection of spaces and mappings as follows; 

K" +— !- I"xG G 

U (17) 
L 

where R" is a real Euclidean space, G is some locally compact Abelian group, and L C K" x G is a 
lattice, i.e. a discrete subgroup for which the quotient group (M™ x G)/Z is compact. Furthermore, 
we assume that tt\\i is injective and 7r 2 (L) is dense in G. 

A model set in R" is a subset of M™ which, up to translation, is of the form A(V) = { n\(x) | x € 
£,7T2(x) G T^} for some cut and project scheme as above, where V C G has non-empty interior 
and compact closure (relatively compact). When we need to be more precise we explicitly mention 
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the cut and project scheme from which a model set arises. This is quite important in some of the 
theorems below. Model sets are always Delone subsets of R", that is to say, they are relatively dense 
and uniformly discrete. 

o 

We call the model set A(V) regular if the boundary dV — V\ V of V is of (Haar) measure 0. We 
will also find it convenient to consider certain degenerate types of model sets. A weak model set is a 
set in R n of the form A(V) where we assume only that V is relatively compact, but not that it has 
a non-empty interior. When V has no interior, A(V) is not necessarily relatively dense in l n but 
regularity still means that the boundary of V is of measure 0. 

Theorem 2 (Schlottmann Q/J If A = A(V) is a regular model set, then A is a pure point diffractive 
set, i. e. the Fourier transform of its volume averaged autocorrelation measure is a pure point measure. 

□ 

It is this theorem that is a prime motivation for finding criteria for sets to be model sets. 
Now let (U, $) be a substitution system with inflation Q on a lattice L of R n and let L be the 
Q-adic completion of L. This gives rise to the cut and project scheme. 





R" x L - 


^ I 




u 




L < — 


L 


-> L 




(*.*) - 


-» t 



(18) 



where L := { (t, t) \ t e L] C R" x L. 

We claim that (W 1 x L)/L is compact. L is clearly discrete and closed inR"xL. Since (R n xL)/L 
is Hausdorff and satisfies the first axiom of countability, it is enough to show that it is sequentially 
compact ||. If {(xi,Zi) + L} is a countable sequence in (R n x L)/L, then there is a subsequence 
{(xi,Zi) + L}s with {xi + L}s convergent sequence, since W 1 / L is compact. We can rewrite 
{(xi,Zi) + L}s as {(x-,z,-) + L}s, where {x'^i^s converges to x in R™. Since L is compact, there 
is a convergent subsequence {^}s' to some z in L. Thus {(x^, z[)} s' converges to (x,z) in R™ x L. 
Therefore (R n x L)/L is sequentially compact. 

Note also that 7Tx | ? is injective and TT2(L) is dense in L. 

Lemma 5 Let Ui, i = 1, — , m, be disjoint point sets of the lattice L in R™. Identify L and its image 
in L. Let Wi := Ui in L and Ai := Wi PI L. Suppose that fi(dWi) = for all i = 1, • • • , to. 

. o 

(i) If Ai\Ui C Uj=i 9Wj then, relative to the CPS (jlq,), C/i is a regular weak model set when Wi 

o 

is empty, and Ui is a regular model set when Wi is non-empty. 

(ii) If L = [Jj—i Uj and each Ui is a regular model set, then Ai\Ui C Uj=i &Wj for all i = 
1, . . . ,m. 

PROOF: (i) Assume that A t \U t C \J™ =1 dWj for all i = l,---,m. Since /j,(dWi) = for all i = 
1, • • • ,TO, 

m 

fi{W t ) = n{Wi) = \ |J dWj) (19) 
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Since A* = Wi H L, Ui = Vi n L where Vi := W l \(A l \U t ). Now V t DW t \{JT=i dW j- From 
Wi \ Ujli C ViC Vi C W = Wi and (|l|), /i(VV\ = 0. So £/, is regular. If Wi= 0, then V= 

o o 

also. Thus Ui is a regular weak model set. On the other hand, for any i with 0, and Vi 
is compact. It follows that Ui = A(V^) is a regular model set for the CPS (|l8|). 

(ii) Suppose that Vrf 0, fi(V t \ V) = 0, where Ui = Vi H L, and i = Ujli tj- Then from 

AAC^ = A(Wi)\A(Vi) C C Wi\ W\ A. we have //(A~\E/j) = for alii = 1, • • • , to. By 

Lemma | (i) and (ii) , Wi n for all t / j and A^ C UJli □ 

Theorem 3 Lei (J7, $) oe a primitive substitution system with inflation Q on the lattice L in K™. 
Suppose that PF-eigenvalue of the substitution matrix S{<&) is equal to |detQ| and L = U i=1 Ui- 
Then the following are equivalent. 

(i) There is a primitive substitution matrix '5 admitting a coincidence, where (U, vp) is equivalent 
to (U, $ M ) for some M > 1. 

(ii) The sets Ui, i = 1, • • • , to, of f/ are model sets for the CPS fildj). 
('wj For ai /easi one i, £/j contains a co set a + Q M L. 

(iv) (U,$>) admits a modular coincidence. 

PROOF : 

(i) =>• (ii): Suppose that (U, \l/) admits a coincidence and is equivalent to (U, $ M ). Fix i £ 
{1, • • ■ ,to} with H^Lj^jj ^ and let g be in this intersection. Recalling equation (^), and in view 
of the choice of g, we have 

n(Wi)< (jr J] a*(/(W,0)) -^( 5 (^)n 3 (^)), 

for any fc , / G {1, • • • , to} with k ^ I. On the other hand, from Theorem [l] (ii) 

m m 

A»(Wi) = — EW^/iTO E M(/(^-)) ■ (20) 

o o 

Thus, in fact, /i(g(Wfc) n g{Wi)) = whenever k ^ I. It follows at once that Wk f~l W;= for all 
k ^ I, since the measure of any open set is larger than 0. 

Recall that Wij^ and fx{dWi) = for all i = 1, • • • , to . Then by Lemma |J(ii) and Lemma 
t/j, i — 1, • • • , to, are model sets in CPS(|l8|). 

(ii) (iii): Assume that Ui, i = 1, • • • , to, are model sets in CPS(|T^), i.e. Ui — A(T^) — Vi D L 

o o 

for some Vi with 0. Thus there is a coset a + Q M L cVi and, since we can always choose the 
coset representative from the dense lattice L, we can arrange that a + Q M L C Ui. 

(iii) =>• (iv): Assume that for at least one i, Ui contains a coset a + Q M L. Fix i. Iterate 3> 
M-times. Then each function / in the substitution system <E> M has the form / : x i— > Q M a; + b. For 
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each j, let Gj := { f G (& M )ij \ t(f) = a mod Q M L}. (Recall that t(f) is the translational part of 
/). From Ui = U™ i \J fe( ** )tj f<Pi), we obtain a + Q M L C |J™ i U /eGj f(Pj). In fact 

m 

a + M L=|J |J /([//), (21) 

since the right hand side is clearly inside a + Q L. From the fact a + Q AI L C Ui, we get $ M [a] = 
Uj=i C Uj=i )w ■ Therefore $ M has a row containing an entire congruence class $ M [a]. 

(iv) => (i): Assume $ M has a row, say i-th row, containing an entire congruence class $ A/ [a]. 
Let Gj := $ M [a] n ($*%•. Then (J£Li U/ec, /(^j) c a + <9 Mi - Reca11 that Ujli ^ = L and 
17 = §> M {U). It follows that the elements of a + Q M L can be obtained from the substitution system 
<1> M only from the mappings of <1> M [a] , and indeed they must all appear as images of the mappings 
of $ M [a]. Thus 

m 

a + Q M L=\J |J f{U s )<zUi. (22) 
3=1 feGj 

On the other hand, 

m 

a + Q M L=[jQ M (U j ) + a, (23) 
j=i 

which is a disjoint union. 

We now alter our substitution system <I> M as follows: Define g : L — > L by g(x) = Q M a; + a. We 
may, by restriction of domain, consider g as a function on JTj, j = 1, • • • ,m. We define 4" by 

f * y = (($ M )^ ag.ouw 
1 * fej = ($ M ) fej if fc^t, 

for all j. From ( |22] ) and (|23]), the ^y, j = 1, • • • , m, consist of maps from Uj to Ui and have the 
same total effect on Ui as the (<& M )ij ,j = 1, • • • , m . Thus (U, ^f) is a substitution system admitting 
a coincidence. 

Since 5($ M ) is primitive, the incidence matrix /(<!> M ) is primitive. Then I(^>) is also primitive, 
since 7(<f> M ) < So $ is primitive. In addition, ^ has the inflation Q M for L which is an 

inflation in $ M . 

We claim that S^), S($ M ) have the same PF-eigenvalue and right PF-eigenvector. Then (U, 4') 
is equivalent to (U, $ M ). 

o o 

We verify first that Wk f~l W}= for all fc 7^ j. We can assume that m > 1, since there is 
nothing to prove when m = 1. Let .91 € Gi = ($ M )iz[a] ^ for some Z. Take any /c € {1, • • • , m}. 
There is M £ Z + for which ($ M °)» 7^ 0. Choose / <E ($ M °)ik- Let g x : x i-> <3 M £ + ai, where 
01 = a mod Q M L, and / : x i-> Q m °x + 6 with tel. Then g x of ix\-+ Q M+M °x + Q M b + at. So 
51 o / G ($ M+Mo ) lfc [ai + Q M b}. Furthermore (o x + Q M 6) + Q A/+M °(L) Ca t + Q M i C C/ 4 . 

Let TV := M + M , c := a a + Q A/ &, andp := 51 o /. Note that 

m 

c + Q N L=\J |J M^), (24) 
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where Hj = ($ w ) y [c]. 

There are at least two functions in Ujli since for all j Uj, ^ L. We can write c + Q N L in 
the form 

c + Q N L = \J{ Q N U 3 + Q N a h + c | j e {1, • • • ,m}, h G Hj,a h G L } , (25) 

where we have used the explicit form of each of the mappings h 6 Hj. This union is disjoint, and as 
a consequence the elements ah £ L for h in any single Hj are all distinct. In particular we have a p 
coming from H^. From (|2^) we have 

m 

i=U U (Vi+oth) (26) 

3=1 heHj 

and separating off 

m 

L = U k u\J {J{Uj+a h - otp) , (27) 

3=1 fcefl* 

where i?j := Hj if j ^ k and -ff^, := Hk\{p}. Again these decompositions are disjoint. But we also 
know that Uk and |Jj=i are disjoint, and it follows that 

m m 

\JUjC\J |J ah-op). 

3=1 j=lhEH' 

Taking closures, 

m m 

U^ C U U (Wj+a h -a p ). (28) 
On the other hand, if we apply Theorem [l](ii) to $ N and look at ( p4[ ) we see that 

m m 

M ( c + q^l) =2 E "WO) = E E MQ^- + + c), 

3=1 heHj 3=1 /iSiJj 



and hence 



Thus 



M L ) = E E ^(^i + ah ) = E E K^i + - Op)- 

3=1 heHj 3=1 ftGi/j 



= M (w fe ) + ( E E /*w + ^ - a p) 

i j=i heH' 



which, after taking closures in d27n, gives 



/» ( ^ n ( U U w> ■ + ^ - a p) ) I = • ( 29 ) 

i j=i hem 
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Finally from (28) and (£9|) we obtain 



n(w k n({JWj)) = o, 



3 = 1 



from which Wt H Wj= for all k 7^ j. This establishes the claim. 
Now 



from = 0, Win Wj= for all j ^ j 

EmsTO). ( 3 °) 



Again using Theorem [l] (ii), this time for <£> M , we obtain 

1 .SY-1' U : 



|detQ M | 

where w = \fi(Wi), . . . , [i(W m )] T . The part of this relation in Wi which pertains to the coset a+Q AI L 
is 

m 

fi(a + Q M L) = E E TO) ■ ( 31 ) 

But from (||) 

Ai(a + Q M L)= M ( \)g{W )\ • (32) 



Together, (RG), fey), and (B2J) show 



1 

IdetQ^ 



Since w > and S(^>) is primitive, S^W) has PF-eigenvalue |detQ M | and PF-eigenvector w as 
required. □ 

Remark: Let A — {ai, . . . ,a m } be an alphabet of m symbols and let cr be a primitive equal- 
length alphabetic substitution system on A, that is, 

(i) a : A — > A q for some q £ Z + ; 
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(ii) the m x m matrix S = (SV/), whose i,j entry is the number of appearances of in 17(0^), is 
primitive. 

According to Gottschalk R], for some iteration a k of a, there is a word w € A z which is fixed by 
er in the sense that 

a k (w Q Wi . . .) = io tWi... , 33 > 
er fc (. ..io_ 2 w_i) = . ..tu_ 2 iw-i • 

Replacing cr' c by <t and q k by g if necessary we can suppose that k = 1, and assume then that 
cr(u>) = w. 

We can view w as a tiling of R by tiles of types 01, . . . a m , all of the same length 1. If we 
coordinatize each tile by its lefthand end point so that wi gets coordinate I, then we obtain a 
partition U\ U . . . U U m of Z and an m x m matrix substitution system $ of g-affine mappings derived 
directly from a: namely, aaj = .. .ai gives rise to the mappings (x 1— > qx + I — 1) G ^.j) 
/ = 1, . . . ,q. 

We take as our cut and project scheme 



R <- 


- RxZ, - 


-» z 9 




u 




Z <- 


z 


z 


z < — 


- ( z , z ) - 


— » z 



(34) 



(see |H), where Z g is the g-adic completion of Z. 

According to Theorem ||, the t/i are model sets for ([m]) if and only if for some iteration a M of a, 
there is a fc G Z for which all the mappings fi : x 1— > q M a; + ^ with ^ = fc (mod g M ) lie in one row of 

Since a M a,j has q M letters in it, there are q M mappings in the jth column of 4> M . Furthermore, 
since the letters a M a.j are represented by contiguous tiles, their coordinates fall in a range of consec- 
utive integers, and so the mappings of the jth column of $ M are the maps /; , where < I < q M , in 
some order. In particular, all of the mappings in $ M are of this restricted form. It follows that mod- 
ular coincidence is equivalent to the existence of a row of $ M , say the ith row, and a fc, < fc < q M , 
so that fk belongs to each of <&,f/, . . . , 

This condition precisely says that there is a fc so that the fcth position of a M (a,j) contains the same 
letter a, for all j. This is the well-known coincidence condition of Dekking ||, and he has proved 
that for non-periodic primitive equal-length substitutions, this condition is equivalent to pure point 
diffractivity. It is straighforward to show that S'( < ^ , ) has its PF-eigenvalue equal to | detQ|. Thus we 
have 

Corollary 1 Let a be a primitive equal-length (— q) alphabetic substitution with a fixed bi-infinite 
word w, and assume that w is not periodic. Let $ be the corresponding matrix substitution system 
and let Z = U± U . . . U U m be the corresponding partition of Z. Then the following are equivalent: 

(i) there is an M so that o~ M has a coincidence in the sense of Dekking; 

(ii) $ has a modular coincidence; 
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(Hi) the Ui 's are model sets for (pjj); 
(iv) the Ui 's are pure point diffractive. 



We note that this interesting equivalence of model sets and pure point diffractivity is more than 
we can yet prove in the higher dimensional substitution systems. 

5 Sphinx tiling 

Long we sought the wayward lynx 
And bowed before the subtle sphinx 
But solved we not the cryptic sphinx 
Before we found the wayward links. 

-Anon 

In this section we take up the sphinx tiling. This is a substitution tiling whose subdivision rule 
is shown in Figure and Fi gure 4 




Figure 1: Sphinx Inflation [Type 1] 




Figure 2: Sphinx Inflation [Type 2] 

It has 12 sphinx-like tiles (up to translation). If we choose a single point in the same way in each 
sphinx then we arrive at 12 sets of points. We wish to show that each of these sets is a regular model 
set. Actually we make a slight alteration to this, choosing several points from each tile, but this is 
equivalent to our original problem. 

Each sphinx can be viewed as consisting of 6 equilateral triangles of two orientations. In this 
way, any sphinx tiling determines a tessellation of the plane by equilateral triangles. We consider 
the centre points of the triangles of one orientation. These clearly form a lattice L, once we have 
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chosen one of them as the origin. Note that some sphinxes have two points and others have four 
points in L. We give names to each tile and the points in it as shown in Figure]^. Then the 12 types 
of sphinx partition L into 36 subsets forming a matrix substitution system. We show that these are 
model sets for a 2-adic-like cut and project scheme of the form of (|Tg|). 
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Figure 3: 12 Sphinx Tiles 

With the origin as shown, the coordinates are chosen so that in the standard rectangular system 
(1,0) is the lattice point directly to the right of (0,0). It is more convenient to replace this by an 
oblique coordinate system: L — { ae + bw \ a, b G Z}, where e = (1, 0), w — (4, ^) in the standard 
rectangular system and relative to this basis we can identify L and Z 2 and denote ae + bw by (a, b). 
The basic inflation shown in Figure [l] gives rise to the map 

T : x i-> 2Rx+ (1,0) , 

where R is a reflection in R 2 through x-axis, i.e. in the new coordinates, R(l, 0) = (1, 0), R{0, 1) = 
(1,-1). 

The various types of points are designated by letter pairs ia, where i £ {1, • • • , 12} and a £ 
{a, ■ ■ ■ , d} (of which only 36 actually occur). Let Ui a be the set of points of type ia. On the basis 
of this we can make mappings of each point set to other point set. 
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Define 

hi : x i-> Tx + (0,0), /i 2 : a; h-> Tx + (1,0) 

/13 : a; 1— > Tx + (0, 1), /14 : x 1— > Tx + (—1, 1) 

: x 1— ► Tx + (—1, 0), h§ : a; 1— * Tx + (0, —1) 

h 7 : x i-*-Tx+ (1,-1), hg : x Tx + (2, — 1) 

Tig : x 1— > Tx + (—1, 2), /iio : x 1— > Tx + (—1, — 1). 

Let fi a jp be the function which maps j/3-point set into ia-point set. 
[Type 1] 



,/*9a la. 


/14 ; 


x 1 — 


Tx + 


f-1 1) 


fla lb 




X I — 


Tx + 


fl 0) 


f$b la 


/lQ • 
f . 


X 1 — > 


Tx + 


f-1 2) 


fib lb 




X I — > 


Tx + 


fo 0) 


/9c la = 


^3 : 


x 1— » 


Tx + 


(0,1), 


/4a lb — 


h 5 : 


X I — 


Tx + 


(-1,0) 


fvd la — 


h 2 : 


X 1 — 


Tx + 


(1,0), 


fiblb — 


hio 


X H 


-> TzH 


-(-1,- 


fia la 


fcl : 


X 1 


Tx + 


(0,0), 


ficlb — 


he ■ 


X I— » 


Tx + 


(0,-1) 


fib la = 


/i 6 : 


X 1 * 


Tx + 


(0,-1), 


fidlb — 


h 7 : 


X I— 


Tx + 


(1,-1) 


/4c la — 


h 7 : 


X 1 > 


Tx + 


(1,-1), 












Hdla — 


h 8 : 


X 1 — 


Tx + 


(2,-1). 












J 12a 4a 


= Ai 


: x 


1— > Tx 


+ (0,0), 


/la 46 


= h 2 


: x 


— > Tx 


+ (1,0) 


f 12b 4a 


= hi 


: x 


— * Tx 


f (-1,1), 


fib 4b 


= hi 


: x 


— > Tx 


+ (0,0) 


fiaic - 


--hi 


: x h 


-> Tx4 


-(0,0), 


fla4d 


= hi 


: x 


i— > Tx 


+ (0,0) 


fib Ac ~ 


-■ h 6 


X H 


•> Tx + 


(0,-1), 


flb4d 


= h 5 


: x 


— » Tx 


+ (-1,0 


fie Ac — 


-h 7 


X 1- 


4 Tx + 


(1,-1), 












fid 4c = 


= h 8 


X H 


-> Tx-f 


(2,-1). 













1) 



All points in a sphinx having 2-points in it are mapped as in [Type 1] changing the translation 
part according to the orientation of the sphinx relative to sphinx 1. Likewise, all points in a sphinx 
having 4-points in it are mapped as in [Type 2] relative to sphinx 4. 

Now we can list the 36 x 36 matrix($) of affine mappings that make up our substitution system 
(Figure |). 

We can check that S , ($) has PF-eigenvalue 4 and is a primitive matrix and the union of point 
sets is L. We used Mathematica to check that property (iv) in Theorem || is satisfied in $ 8 (it may 
actually be satisfied at some lower power). Certainly in <I> 8 there are a large number of modular 
coincidences. Theorem [l] and || say that all 36 point sets are regular model sets in CPS (|l8|). 

6 The total index and model sets 

In this section we derive another criterion for determining when a partition of a lattice is a partition 
into Q-adic model sets, the difference this time being that there is no substitution system involved. 
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((0,h2,0,0,0,0,0,0,0,h2,0,hl,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,hl,0,0,0,0,0), 
(0, hi, 0, 0, 0, 0, 0, 0, 0, hi, 0,h5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h5, 0,0,0,0,0}, 
{0,0,0,0,0,0,hl,h7,0,0,0,0,0,0,0,0,hl,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h7,0), 
{0,0,0,0,0,0,h2,h8,0,0,0,0,0,0,0,0,h2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h8,0), 
(0,0,0,0,0,0,h3,h2,0,0,0,0,0,0,0,0,h3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h2,0), 
(0,0,0,0,0,0,h9,h3,0,0,0,0,0,0,0,0,h9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h3,0), 
(0,0,0,h6,0,hl,0,0,0,0,0,0,0,h6,0,0,0,0,hl,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0), 
(0,0,0,hl,0,h3,0,0,0,0,0,0,0,hl,0,0,0,0,h3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0), 
(hl,hS, 0,0, 0,0, 0,0, 0,0, hi, 0,0, 0,0, 0,0,0,0,0,0,0,115, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0), 
(h6,hl0,0,0,0,0,0,0,0,0,h6,0,0,0,0,0,0,0,0,0,0,0,hl0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0), 
(h7,h6,0,0,0,0,0,0,0,0,h7,0,0,0,0,0,0,0,0,0,0,0,h6, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0), 
(h8,h7,0,0,0,0,0,0,0,0,h8,0,0,0,0,0,0,0,0,0,0,0,h7, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0), 
(0,0,0,0,0,0,0,h4,0,0,0,0,0,0,0,h4,0,hl,0,0,0,0,0,0,hl, 0,0, 0,0, 0,0, 0,0, 0,0,0), 
(0,0,0,0,0,0,0,hl,0,0,0,0,0,0,0,hl,0,h7,0,0,0,0,0,0,h7,0,0,0,0,0,0,0,0,0,0,0), 
(0, 0, 0, 0, hi, 0,0,0,0,0,0,0,hl,h3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h3, 0,0, 0,0, 0,0,0), 
{0,0,0,0,h4,0,0,0,0,0,0,0,h4,h9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h9, 0,0, 0,0, 0,0,0), 
(0,0,0,0,h5,0,0,0,0,0,0,0,h5,h4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h4, 0,0, 0,0, 0,0,0), 
(0,0,0,0,hl0,0,0,0,0,0,0,0,hl0,h5,0,0,0,0,0,0,0,0,0,0,0,0,0,h5, 0,0, 0,0, 0,0,0), 
(0,0,0,0,0,0,0,0,0,0,0,0,h2,0,0,0,0,0,0,0,hl,0,0,0,0,0,0,0,hl,h2,0,0,0,0,0,0), 
{0,0,0,0,0,0,0,0,0,0,0,0,h8,0,0,0,0,0,0,0,h7,0,0,0,0,0,0,0,h7,h8,0,0,0,0,0,0), 
(0,0,0,0,0,0,0,0,0,0,0,0,h7,0,0,0,0,0,0,0,h6,0,0,0,0,0,0,0,h6,h7,0,0,0,0,0,0), 
(0,0,0,0,0,0,0,0,0,0,0,0,h6,0,0,0,0,0,0,0,hl0,0,0,0,0,0,0,0,hl0,h6,0,0,0,0,0,0), 
(0,0,0,0,0,0,0,0,0,0,0,0,0,0,hl,0,0,0,0,0,0,0,0,0,0,h3,0,hl,0,0,0,0,0,0,0,h3), 
(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, h6, 0, 0, 0, 0, 0, 0,0, 0,0,0, hi, 0,h6, 0,0, 0,0, 0,0,0, hi), 
(h4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0, 0,0,0, hi, h4, 0,0, 0,0, 0,0, 0,0, hi, 0,0,0), 
{h9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h3,h9,0,0,0,0,0,0,0,0,h3, 0,0,0), 
{h3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h2,h3,0,0,0,0,0,0,0,0,h2, 0,0,0), 
(h2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h8,h2,0,0,0,0,0,0,0,0,h8, 0,0,0), 
(0,0,hl,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h5,0,hl,0,0,0,0,0,0,0,h5,0,0,0,0,0,0), 
(0, 0, h2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0,0, hi, 0,h2, 0,0, 0,0, 0,0,0, hi, 0,0, 0,0, 0,0), 
(0, 0, 0, 0, 0, 0,ht, 0, 0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, hi, 0,0,0,0,0,0,0,U.,M), 
(0,0,0,0,0,0,hl0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h5,0,0,0,0,0,0,0,h5,hl0), 
(0,0,0,0,0,0,h5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h4,0,0,0,0,0,0,0,h4,h5), 
(0,0,0,0,0,0,h4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h9,0,0,0,0,0,0,0,h9,h4), 
{0,0,0,0,0,0,0,0,hl,0,0,0,0,0,0,0,0,0,0,0,0,0,0,h7,0,0,0,0,0,0,0,h7,0,hl,0,0), 
(0, 0, 0, 0, 0, 0, 0, 0, h4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0, hi, 0,0, 0,0, 0,0,0, hi, 0,h4, 0,0)) 



Figure 4: Sphinx matrix function system(<f>) 

We assume that we are given a lattice L in R" and an inflation Q on L as in (||). The notation 
remains the same as before. The main ingredient is a non-negative sub-additive function called the 
total index which is defined on the subsets of L and its Q-adic completion L. 

For any subset V of L the coset part of V is defined as 

C(V):=|J{C|Cisacosetin^}. (35) 

The key point to remember in what follows is that two cosets in L (L) are either disjoint or one of 
them is contained in the other. If C = a + Q k L is a coset then we write [L : C] for the index of the 
subgroup Q k L in L. 

Lemma 6 The coset part ofV can be written as a disjoint union of cosets in V. 

PROOF: If V contains no cosets, then the result is clear. Suppose V contains cosets. Let C\ = 
»i + Q kl L be a coset in V with k\ minimal. Consider V\C\. No coset can be partly in C\ and 
partly in V\C\. Thus, if V\C± contains no cosets, then C(V) = C\. Otherwise let C2 be a coset 

a 2 + Q k2 L with &2 minimal in V\C\. Then C(V) D C\ (J C%. We continue this process. Since 

there are only finitely many cosets for Q k L in L, either we obtain C(V) = C\ [j ■ ■ ■ (J C r for some 

r or C(V) D C\ [j C% (J • • • , where k\ < k 2 < • • • is infinite and unbounded. In the latter case, 
Cty) = Ufci is our required decompostion. If not, there is a coset C = a + Q k L in V such that 
C (t Ui^i Cj- Then there is Cj with < k < fc,. This contradicts the choice of Cj. □ 
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For V C L, we call a decomposition C(V) = (J- C, of C(V) into mutually disjoint cosets using 
the algorithm of Lemma ^, an efficient decomposition of V into cosets. In this case we call c(V) :— 
J^JL : Ci}^ 1 the total index of V. Since any coset is an efficient decomposition of itself, we have 
C (X) = Si c (Ci)- We will see shortly that the total index is finite. 

It is useful to note that an efficient decomposition of C(V) = (J^ Ci of C(V) into cosets has the 
following special property: if D is any coset of V then necessarily D C Ci for some i. 

Lemma 7 Any two efficient decompositions of C(V) are the same up to rearrangement of the order 
of the cosets. In particular the total index is well-defined. 

PROOF: Let C(V) = [j C be a second decomposition of C(V) determined by the same algorithm as in 
Lemma ||. Then with k\ as in the Lemma, let D%, ■ ■ ■ , D r be all the cosets of V of the form a + Q kl L. 
These are all disjoint and by the algorithm all of them must be chosen in the decomposition of C(V), 
and they all occur before all the others. Thus C\, ■ ■ ■ , C r and C[, ■ ■ ■ , C' r are Di, • • • , D r in some 
order. Removing these and continuing in the same way the result is clear. □ 

We have similar concepts in L. For W C L we have the coset part C*(W) of W and C*(W) can be 
written as a disjoint union of cosets in W. Let C*(W) = lj i Di where Di, i = 1,2, • • • , are mutually 
disjoint cosets in W . We call c*(W) :— ^2 t [L : Di]^ 1 the total index of W. This time we do not need 
to be careful about the way in which the decomposition is obtained since the total index is nothing 
else than the measure /x(C*(W)) of C*(W). 

Given an efficient decomposition C(V) = Ui=i Q into disjont cosets in L, we define C(V) := 
Uj=i Ci C L. This is actually an open set in L. Since [L : C] = [L : C] we see that c(V) = c*(C(V)). 
In particular it follows that the total index of any subset V of L is finite and bounded by n(C(V)). 

Lemma 8 For X,Y C L and X C Y, and any decomposition C{X) = |J i Ci into disjoint cosets, 
E 4 c (C0 < c(Y). In particular, c(X) < c(Y). 

PROOF: Assume first that Y is a single coset C. Then 

]T c{C t ) =Y J c* (Ci) - ]T pQOi) < li{C) = c* (C) - c(C) , (36) 

i i i 

since the cosets remain distinct after closing them in L. 

In the general case, let C(Y) — |Jj=iCj be an efficient decomposition of Y. Since X C Y, each 
Ci C Y. In view of the remark above about efficient decompositions, there is for each i a unique j 
for which d C Cj . Thus we can arrange the Ci 's so that 

oo 

C(X) = |J |J a (37) 

where Aj := {« | C; C Cj}. NowlJ ieA . Cj C Cj, so by the first part of the proof, X^eA c(Ci) < c(Cj). 
Finally 

< x ) = E E c (^) ^ E c (^) = c ( r ) ■ ( 38 ) 

□ 
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Lemma 9 Let Ui, i = 1, ■ • • ,m, be disjoint point sets of the lattice L in W 1 . Let Ai = Ui fl L and 
C(Ui) be the coset part in Ui. Then Ui=i(^A^) C L\\J r ?L 1 C(Ui) ) with equality if L = (J™ x Ui. 

PROOF For x G lj"=i C{Ui) there is a coset C C C(t/j) for which ieCc [/*. Let C = a + Q k L, a e L. 
Suppose a; is a limit point of Uj in L for some j 7^ i. Then, since a + Q k L is an open neighborhood of 
x, (a + Q k L) H Uj 5^ i.e. (a + Q fc L) n J7, 7^ 0. But then U t n f/j 7^ 0, contrary to the assumption. 
This means x g" Ul=i(^A^)i proving the first part. 

Suppose that L = (J™ 1 f/j and x £ L but x ^ U"=i C{Ui). Then x £ Ui for some ?7j but there is 
no coset in [7, which contains x. For any k G Z+, Bk(x) := x + Q k L is an open neighborhood of x 
in L and L n Bk{x) <£_ Ui, by assumption. Since L — \J^ 1 Ui, (L n Bk(xj) fl Z7,- 7^ for some j 7^ i. 
So we can choose x\ G (L n -Bfe(x)) (~l E/j. Then we get a sequence {x£} convergent to x as k — > 00. 
Choosing a subsequence lying entirely in one Aj shows that x G Aj for some j 7^ i. Since x G t/j, 
and Ui,Uj are disjoint, x G Aj\Uj. □ 

Theorem 4 LetUi, i — 1, ••■ , to, be disjoint nonempty point sets of the lattice L inW 1 . LetC(Ui) be 
the coset part in Ui, c(L^) the total index ofUi, and Wi the closure ofUi in L. Then 53j=i C (£A) = 1 
if and only if the sets Ui, i — 1, • ■ ■ , to, are regular weak model sets in the CPS |lS|) and L = (JjJLi W*- 

PROOF 

(=>) Assume that X)"=i c (^i) = 1- Let t/ m +i := L\\J£ =1 U m - Using Lemma ^ and the fact that 
c(L) = 1, we see that c(U m+ \) = and Y^T=i c (Ui) = L For this reason we can assume, in proving 
that the {/, are weak model sets, that Ul=i Ui = L in the first place. 

For j ^ k the cosets of C(Uj) (of which there may be none!) and those of C(Uk) are disjoint from 
one another, and the same applies to C(Uj) and C{Uu). Thus 



n U = E = E °w = 

\i=l / i=l 

and 



M^iVlJC^)) J =0. (39) 
Now note that dVF, n \J^ =l C{Ui) = for any j. If not let a G nC(£/ fe ) for some k. Since 

o 

C(Uk) CWk, we see that j 7^ fc. But a G W}, so a is a limit point of Uj ,and C(Uk) is an open 
neighborhood of a, so t/j fl C(Uk) 7^ 0. This violates the disjointness of the Ui's. We conclude that 
dWj C I\(U™iC(C/i)) and hence that 

K9Wj) = 0, (40) 

for all j = 1, ... , to. Note also 

m 

AAUi c LlCWi) 

i=i 

m 

C L\ [J C(fZj) by Lemma 8 
j'=i 
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L\\J(c(u j )nL) = L\\Jc(u j ). 

3=1 i=i 



This shows that 



By Lemma | (i) and (ii), W, n Wj= for all i ^ j and (A^i) C Ujli ^Wj for all i = 1, • • • , m. 

Using Lemma ||(i) we obtain that the sets U,i = 1, ■•■ ,m, are regular weak model sets in the 
CPS©. 

o 

Remark : Whenever Wi^ 0, U is actually a regular model set. 

Since \J" =1 C{Ui) c U"Li W *> KU?=i W *) = 1 - Thus AU™i ff - is °P en of measure and 
L = UI=i Wi- This last argument does not require that U"=i U = L. 

o 

(<=) Assume that Ui = A(Vi) — Vi fl L where Vi\ Vi has measure and L = {J™=iWi. Thus 

o 

Ui C Vi and Wi := Ui C Vj. Since L is dense in L and for x QVi each ball around a; of radius 

o o o 

e > contains points of Vi f\L C C/j, it follows that Ui DV*. This proves that ViC Wi C Vi. So 
= n{Wi) and = 0. Now 

m rn m 

{JW i = ({Jv i )U(\J(Wi\V i )). 

i—1 i—1 i—1 

o o o 

So ^(UHi Wi) — m(U£Li K)« Also the disjointness of the Ui gives D Vjf= for i ^ j (since L is 
dense in L). Finally 

m m m m m 

2 — 1 Z— 1 Z— 1 i — 1 Z— 1 

□ 

Corollary 2 Let (J7, <&) be a primitive substitution system with inflation Q on the lattice L in R n . 
Suppose that PF-eigenvalue of the substitution matrix S($) is equal to |detQ| and L = (Jt=i U- 
Then Y^Li c(Ui) = 1, where c(U) is the total index of U if and only if the sets Ui, i = 1, • • • , m, are 
model sets in CPS (fT8|) . 

o 

PROOF Use Theorem Ffl to determine that for all i, Wi^ 0. Now use Theorem and Remark in the 
proof. □ 



7 Chair tiling 

The two dimensional chair tiling is generated by the inflation rule shown in Figure |^. There are 
4 orientations of the chairs in any chair tiling. In |0] it was shown that the chair tiling has an 
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interpretation in terms of model sets based on the lattice Z 2 and its 2-adic completion as internal 
space. 

In this section we generalize this result to the n-dimensional chair tiling using the results of the 
last section (see Figure for an example of the 3-dimensional chair). To make things clearer we 
begin with the case n = 2. 



Figure 5: 2-dimensional chair tiling inflation 



I. Chair tiling in R 2 

The starting point is to replace each tile by 3 oriented squares. Figure ^ shows the inflation rule, for 
one chair, in terms of oriented squares. The resulting tiling is a square tiling of the plane in which 
each of the squares has one of 4 orientations. The centre points of each square form a square lattice 
which we identify with Z 2 by assigning coordinates as shown. 
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Figure 6: 2-dimensional chair tiling substitution 



Let Ui be the set of centre points corresponding to squares of orientation (i) as given in Figure |^. 
We start out from a basic generating set A 2 := { (xi,X2) | Xi G {0, —1} } and determine the precise 
maps for the substitution rules of Figure ^. 

Letting e\ := (0, 0), e 2 ■= (1, 0), e 3 := (1, 1), e 4 := (0, 1) , these maps are defined as: 

fj,i ■ Ui -y Uj by (x, i) ^ (2x + ej,j) if j i ± 2 
fiA-Ui^Ui by (x,i) ~ (2x + e,,i) if j = i ± 2 , 
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where i, j e {1,2, 3, 4}, a; e Z 2 , i±2:= 



i + 2 if i < 2 
i - 2 if i > 2 . 



These are the maps of an afhne substitution system $. In fact, if we define 

h\ : x i— > 2x + ei, /12 : x i— > 2a; + e 2 , /13 : x h- > 2.x + e 3 , /14 : x 2a; + 64, 

then 

f\,i = h\, fi t 2 = hi, f[^l = h 3 , fi,4 = hi 

h,l — hi, f 2 ,2 = /2,3 = ^2, = ^4 

/3 2 3 = h l, /3,2 = ^3, /3,3 = ^3, /3,4 = ^3 
/4,1 = /l4, /ifi = h 2 , fi,3 = hi, / 4j 4 = hi , 

( {hi,h 3 } {hi} {} {hi} \ 

{M {h 2 ,h 4 } {M {} 

{} {h 3 } {h 3 ,hl} {h 3 } 

\ {hi} {} {hi} {hi,h 2 } J 

Inflating A 2 by the substitutions above we generate the 4 point sets Ui, i = 1,2, 3, 4. The precise 
description of Ui is the following : 



and 



(7 2 



U 3 



Ui 



Q (J ((0,0) + 2 fe (2,0)+t(l,l) + 2 fe -4Z 2 ) 
fc=o t=o 

oo 2 fc -l oo 

U (J U ((°- °) + 2 *( > 2 ) + !) + ^ ' 4Z ') U U W 1 ' !)> 

fc=0 t=0 t=-oo 
oo 2 fc -l 

(J (J ((-1, 0) + 2 fe (2, 0) + i(-l, 1) + 2 fe ■ 4Z 2 ) 
fe=0 t=o 

oo 2 fc -l oo 

u U U ((-l,0) + 2 fe (0, 2) +t(-l, l) + 2 fe -4Z 2 ) U |J{(0, 1) +t(l, 1)} 

fe=0 t=0 t=0 
oo 2 fc -l 

|J (J ((-l,-l) + 2 fe (2,0)+t(-l,-l) + 2 fe -4Z 2 ) 
fc=o t=o 

oo 2 fc -l 

U |J |J ((-l,-l) + 2 fe (0,2)+t(-l,-l) + 2 fe -4Z 2 ) 

k=0 t=0 
oo 2 fc -l 

|J (J ((0, -1) + 2 fe (2, 0) + t(l, -1) + 2 fe ■ 4Z 2 ) 

oo 2 fc -l oo 

U U U ((0,-l) + 2 fe (0,2)+t(l,-l) + 2 fe .4Z 2 ) U (J{(-l, ) + i(-l,l)}. 



k=o t=o 



t=0 
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Each of these decompositions is basically into cosets, with the exception of three trailing sets in 
types 1,2,4 which we will designate by V\, V2, V\ respectively. 
We can prove the correctness of this as follows: 

Let U[,U2,U^, U' A be the sets on the right hand sides respectively. Note that 

(i) The generating set A 2 is contained in U[ adequately, i.e. 

(0, 0) e U{, (0, -l) e U 2 , (-1, -1) e U[, (-1, 0) e Ui 

(ii) Claim that U[ D (jj =1 C7j, i = 1, 2, 3, 4. 
Check that for any i, 

4 00 2 fc -l 



Wi) C IJ U U ((-e l ) + 2 fc+1 (2(e l -e,)) + 2t(e l±2 - ei ) 

+ 2 fe+1 • 4Z 2 ) U Vi 



j = i k=0 t=o 



c f/i 

4 00 2 fc -l 

fti±2(t^) C IJ U U ((-e l ) + 2 fe+1 (2(e l -e,)) + (2t + l)(e 4±2 -e l ) 

i=i fc=0 t=0 

+ 2 fe+1 • 4z 2 ) u v; 

^(C70 C (-2e ; + e t +4Z 2 )U(-2e i±2 + ei +4Z 2 ) 
C t/-, where I =fi,i± 2, / 6 {1,2,3,4} 



(iii) i = 1,2, 3, 4, are all disjoint. 

Indeed within each Ui, all the cosets and the non-coset part are clearly disjoint. And two 
cosets or non-coset sets chosen from U- and U'j, where j ^ i, i ± 2, cannot intersect, since they 
are different by mod 2. Futhermore two of cosets or non-coset sets chosen from U[ and U[± 2 
cannot intersect either, since for a + 2 k ■ 4Z 2 C U[, b + 2 l ■ 4Z 2 C U' i±2 with k <l, a -b ^ 
mod 2 k ■ 4Z 2 . 

Now since U'^U^U'^U^ are generated from A 2 by $, C/j C U{ for all i = 1,2,3,4. Also from 
U Ui = Z 2 , we get U £/? = Z 2 . Since all U[, i = 1, 2, 3, 4, are disjoint, C/ 4 = f/? for all i = 1,2, 3, 4. 
Finally, for any i = 1, 2, 3, 4, 

c(Ui) > 2 • X) E ( 2 fc. 4 )2 = 2 ' £ i 6 .(2*)2 = 4 • 
Thus X)t=i c ( C/ '0 = L Theorem | shows that U, i = 1,2, 3, 4, are regular model sets. 
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II. Chair tiling in R 



n 



In this section we are going to generalize the foregoing to the n-dimensional chair tilings for all 
n > 2. The n-chair is an n-cube with a corner taken out of it. The inflation rule, which we spell out 
algebraically below, is geometrically the obvious generalization of the 2-dimcnsional case. 



We transform the geometry by replacing each chair by a 2" — 1 oriented cubes, as before, and 
coordinatize the lattice formed by the centres of the cubes, starting from the basic generating set 
A n :— {{x\, ■ ■ ■ ,x n )\xi £ {0, —1}}- There are 2™ orientations of cubes and hence 2" types of points 
(but only 2™ — 1 of these types appear in the starting set A n ). 

For each k > let (3{k) be the binary expansion e D + e\2 + e 2 2 2 + . . . of k, e; £ {0, 1}. We define 
the basic orientation vectors ei, . . . by 



We determine the sets [/,, i = 1, . . . , 2", of all i— type points in Z™ from the points of the basic 
generating set A n , using the inflation rules below. 

The types of the points of A n are as follows: for x = (x\, . . . , x n ) £ A n , 

when x n = —1 ; 

x £ Ui, for which /3(i — I) = (1, . . . , 1) + x, 

when x n = ; 

if x = (0,...,0), x e U± 

otherwise, x £ f7j +2 n-i, for which (3{i — I) = (I, . . . , 1) — ((I, . . . , 1) + x) . 

The idea of considering our vectors in the form (1, . . . , I) + x is to make it easy to compare them 
with the basic orientation vectors. 

This conforms with what happens when n = 2: there are 2™ — 1 types in the basic starting set 
that are in 2™ _1 — I complementary pairs and 1 pair of vectors (0, . . . , 0) and (—1, . . . , —I) of the 
same type, namely of type 1. 




Figure 7: 3-dimensional chair tile 




e„_i) the binary digits of (3(i — 1) if i < 2 
I — e^o"- 1 if i > 2 



in— 1 



in— 1 



Define 



fj,i ■ Ui Uj by (x, i) i > (2x + ej , j) if j^i± 2 
fif -Ui^Ui by (x, i) ~ (2s + ej , i) if j = i± 2 
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i + 2 n - 1 if * < 2" 1 
i~2 n ~ x ifi>2"" 1 . 
Let $ be the matrix function system. Define hi : x i— ► 2x + ei, i £ {!,-••, 2™}. 



where i,j G {1, • • • ,2"}, t ± 2"- 1 

ioi 



$ = 



{/ii,/i 1+2( „-i)} {/ii} 



1 + 2 71 - 
1 

{} 



{M 



V {h 2 A {fl 2 n} ■■■ {} ••• {h 2 „_ 2n -l,h 2 n} J 



Inflating A n by the maps, we get the precise description of U : 

2™ oo 2 fc -l 



Ui 



U U U ((- e <) + 2 "( 2 ^ - e i)) + *fe±2»-i - ei) + 2 k ■ 4Z") |J y 4 



3 = 1 fc = t = 



where 



Ufc-oo{*( e i±2— 1 - ei)} 



if i = 1 



5 ifi^l + 2"- 1 . 

The equalities can be proved in the same way as in the 2-dimensional case. 
Let U'i be the set of the right hand side in ((f2[) . Note that 



(42) 



(i) The generating set A n is contained in U[ adequately i.e. 

e x e U[ if i = 1 

-e i±2 n-i e [/; if i ^ 1, 1 ± 2"- 1 

-e 1±2 »-i eU[ if i = 1 ± 2"- 1 . 

(ii) Claim that U[ D [jf =1 *yt/J,t = 1, • • • , 2". 
Indeed fori e {l,---,2™} 

2" co 2'-l 

M^) C (J U U ((-e,) + 2 fc+1 (2(e 2 -e,)) + 2t( ei±2 „- 

2 fe+1 • 4Z") U Vi 



i - ei 



3=1 fc=0 t=o 



C U' 



h i±2 n-,(Ui) C |J |J |J ({~e l ) + 2 k+1 {2{e l -e )) + {2t + l){e l±2 ^-e l ) 

+ 2 k+1 ■ 4Z") U V$ 



3=1 fc=0 t=o 
j'^i,i±2"- 1 



2G 



hi{U'i) C (-2e ; + ej + 4Z n )U(-2e i±2 „-i + ei + 4Z") 
C U< , where I ^ i, i± 2 n ~ 1 , I e {1, • • • , 2 n } . 

(iii) U[,i= 1, • • • , 2™, are disjoint. 

Indeed all cosets and a non-coset set in each U[ are all disjoint. And two cosets or non-coset 
sets chosen from \J' i and Up where j ^ i,i± 2"~ 1 , cannot intersect, since they are different 
by mod 2. Futhermore two of cosets or non-coset sets chosen from \][ and f/f , 2 „_i cannot 
intersect either, since for a + 2 k ■ 4Z™ C [7?, 6 + 2' • 4Z" C t/- ±2 „-i with fc < Z, a - b ± mod 
2 fc • 4Z n . 

Now since [/,,',« = 1, • ■ • , 2", are generated from A n by [/, C £7/ for all i = 1, • • • , 2". Also from 
Uilx f/i = Z™, (Jill ^ = Z n . Since all V[, i = 1, • • • , 2™, are disjoint, [/, = U[ for all i = 1, • • • , 2". 



For any i = 1, • • • , 2", 

oo 2 fc -l 

TV 



c(C/j) > (2" - 2) • > \ t^t— r- = (2" - 2) 



1 N v-v 2 k 1 



(2*. 4)" v 7 ^ 2 2 " • (2 k ) n 2 r - 



Thus X)i=i c (^i) — 1- Theorem [| shows that [/j, i = 1, ■ ■ • , 2", are regular model sets. 

To get a model set interpretation of the chair tiling itself we proceed as follows. We observe 
that every arrow points to the inner corner of exactly one chair. Let us label each chair by its 
inner corner point which is at the tip of exactly 2" — 1 arrows. These corner points give us 2 n sets 
Xi, . . . , A 2 ™ according to the type, and all lie in the shift L' = (|, . . . , |) + Z™ of our lattice Z n . Let 
fi,i = 1, . . . , 2 n , be (|, . . . , i) — ej respectively. Then Ui +fi is the set of tips of all arrows of type i 

o 

and Ui + fi = V n (Vi + /;) , for some V, C Z£ for which V l compact, and n{dVi) = 0. Now 

Xi = z/ n ( p| (y, + /,•)) 

which is the required regular model set description of Xi, since 

9( n u ^-+&) 

and + /i)) = for all j = 1, . . . , 2™. 

From this result we can show that if we mark each chair with a single point in a consistent way, 
then the set of points obtained from all the chairs of any one type also forms a regular model set, 
and hence a pure point diffractive set. 
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